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Tumour metabolism is an outstanding topic of cancer research, as it determines the growth rate 
and the global activity of tumours. Recently, by combining the diffusion of oxygen, nutrients, and 
metabolites in the extracellular environment, and the internal motions that mix live and dead cells, we 
derived a growth law of solid tumours which is linked to parameters at the cellular leveP. Here we use this 
growth law to obtain a metabolic scaling law for solid tumours, which is obeyed by tumours of different 
histotypes both in vitro and in vivo, and we display its relation with the fractal dimension of the distribution 
of live cells in the tumour mass. The scaling behaviour is related to measurable parameters, with potential 
applications in the clinical practice. 



Recently we derived a growth law for solid tumours\ where growth depends on the distribution of live cells in 
the tumour mass. The model was suggested by the results of computer simulations^, and it has been 
validated using data from in vitro experiments. In this model the inter-vascular micro-regions of solid 
tumours contain both live and dead cells, and the fraction of live cells is exponentially smaller as the distance from 
the sources of nutrients and oxygen increases. Although the decay length 1 is weakly dependent on tumour size, 
we found that an effective, size-independent X works nearly as well, and that X has values mostly in the range 50- 
150 jim. Notably, this is also the distance from the nearest blood vessel where the interstitial p02 assumes hypoxic 
and anoxic values and pH drops to acidic values in the micro-regions of vascularized tumours^ and corresponds 
to the typical thickness of the viable tumour cell layer around blood vessels^. Here we show how the growth law 
can be combined with basic metabolic parameters to yield a seemingly universal metabolic scaling law for 
tumours, that holds both in vitro and in vivo. 

We start from the tumour growth law, that can be expressed as a differential equation for tumour volume, and 
combines proliferation of live cells in the tumour with the gradual shrinking of dead cells ^ 



dV 
dt 



= {oi + d)V{t)F[V{t)]-dV{t), 
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is the fraction of live cells in the tumour, derived from the assumption of exponential decay of the density of live 
cells, so that the total volume of live cells is Va = F(V)V, and where the parameter a defines the individual cell 
proliferation rate, while the parameter 3 is the shrinking rate of dead cells. The tumour volume is proportional to 
x^, where x is some characteristic length of the tumour, i.e., V = Ax^. In the case of spherical tumours A = 4n/3 ~ 
4.2, and x is the tumour radius, but here we consider the possibility of departure from a spherical shape. Still, we 
assume that the tumour keeps roughly the same shape as it grows, so that x can be chosen as a substitute for radius, 
e.g., as the length of a given chord between two fixed, recognizable surface features of the tumour shape. We note 
that A is one of the factors that set the rate of tumour growth, at least for large tumour sizes, since in that case 
F(y)~3M^/^y~^/^, and the growth term in Eq. (1) depends explicitly on A: 



dV 
'dt'' 



^^3W^\a + 3)[V{t)f'-3V(t) 



(3) 



In the context of this model, the nutrient consumption rate fi of the whole tumour is given by the number of live 
cells in the tumour F{V)V/Vc, times the mean consumption rate per cell c 
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where is the mean tumour cell volume, and rj = clv^ is the mean 
nutrient consumption rate per cell and unit volume. We remark that 
Eq. (4) shows that the nutrient consumption rate interpolates 
between a linear behaviour 

ix^nv (5) 

at small tumour size, and a power law with exponent a = 2/3 

li^3kr]A'''V^'^ (6) 

at large tumour size (see Fig. 1). 

At this point we stress that the rate r] = cIVc is cell-type specific, 
while the area parameter A depends on the kind of tumour (and 
therefore on the cell type), on the particular tumour history, and 
finally on the choice of the characteristic length x. These latter two 
dependences should normally be weak, if we assume that x is always 
close to the diameter of a nearly spherical shape, and that the 
development of a given kind of tumour is roughly the same in dif- 
ferent experimental or clinical settings. 

Results 

We have considered the experimental values of the normalized rate 
jLi = lLi/f] of glucose consumption in tumour spheroids obtained from 
different cell lines, as well as in different histotypes of primary and 
metastatic human solid tumours grown in rats^"^^ (see Methods). 
This normalization factors out the specific metabolic consumption 
of different cell lines. From individual fits of each data set we find that 
they are well described by 

fi 3X]^Y 

^^~32^ + (y/A0'/' 



(7) 



where is the cell-line specific consumption rate, and X]^ and Ay. are 
the parameter values obtained for the corresponding /c-th data set. 
The values k]^ and A]^ obtained in these fits are reported in Table 5 in 
Methods, and we notice that Xy^ varies in a very limited range, while 
differs markedly in in vitro and in vivo tumours. The specific 
values Ak were used to define the doubly normalized rate 



3X{V/Ak) 



3Xz 



(8) 



where k is an effective value of the parameter - which is justified 
by the limited variation found in the fits - and z = V/A. Using 

the redefined rate fiP = fi/T]j^Ak, all individual data points 

(^V-^^ /Ak.fif^/rjj^Ak^ {i-th experimental value of the k-th data set), 

from both in vitro and in vivo tumours, follow one common curve 
An ^ shown in Fig. 2, that depends on the effective parameter 

1, which we estimate from the fit. 

While we discuss the whole selection of data in the Methods sec- 
tion, here it is important to note that in the case of human tumours in 
vivo we must cope with the scarcity of available data. Indeed, for each 
data set we need the metabolic rate - i.e., rate of glucose consumption 
- of individual cells, their average volume and the metabolic rate of 
the corresponding tumours, and the combined knowledge of these 
data is quite rare. Therefore we group the few complete data on 
human tumours in vivo in a single data set, and this produces an 
additional variability which is not accounted for by the model. 
However, we find that the rj's for different human tumours all lie 
within a factor 2 of their average (see the Methods section). While 
this adds to the fiuctuations about the scaling curve, we note that 
logio 2 — 0.3, and that the largest fluctuation about the scaling curve 
is about 0.4 - 0.5 as measured on the log scale of Fig. 2, and therefore 
we see that the largest contribution to the observed fluctuations is due 
to this uncertainty on the value of rj. We also remark that the fluctua- 
tions are partly due to our using the same A^ for all tumours of the k- 
th cell line, while, strictly speaking, A is a parameter related to an 
individual tumour. Since the tumours in our dataset follow rather 
neatly the scaling curve in Fig. 2, we argue a posteriori that the 
contribution to fluctuations due to different A's and to different 
tumour environments is less than about 0.2 on the log scale of the 
figure - i.e., a factor smaller than 1.6 on the total tumour metabolic 
rate. 

Even so, data points seem to follow rather nicely the predicted 
model behaviour, and the curve of Fig. 2 provides an indirect con- 
firmation of the validity of the growth law in real tumours. Moreover, 
there seem to be common mechanisms that drive the growth of solid 




Figure 1 | Plot of the normalized nutrient consumption rate ju = ju/tf 
(\im^) vs. tumour volume V (^m^), as given by Eq. (4) (solid line). Here we 
assume a spherical shape and /I = 100 |am, which is close to the values 
found in experimental tumour spheroids \ The dotted line is the 
normalized rate ju from Eq. (5), while the dashed line is ju from Eq. (6). The 
volume range corresponds to a minimum radius of 5 \im (i.e., 
approximately a single cell), up to a maximum of 2000 [xm. The arrow 
marks the volume corresponding to a nearly spherical avascular tumour 
with a diameter of 1 mm. 




• Human tumours in vivo 

• 9L splieroids 

• Rat1-T1 splieroids 

• IVIR1 splieroids 



V/A (jum^) 

Figure 2 | Doubly normalized glucose consumption rate fi^=fi^{z) vs. 
z{z= V/A). Here we take glucose as representative of all nutrients (see 
Methods). The black line is a single fit of Eq. (8) to all data shown in the 
figure: the fit yields 2 = 102 ± 2 |im, and it is compatible with the values 
found in the analysis of growth curves of tumour spheroids \ The dotted 
line is the linear approximation ju^^Zy Eq. (5), at small tumour size, while 
the dashed line is the approximation fi^^3lz^^^, Eq. (6), at large tumour 
size. Data from human tumours (green circles) include breast, uterine and 
ovarian carcinomas, melanomas, thyroid carcinomas, colon and lung 
carcinomas. 
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Figure 3 | Effective fractal dimension of the set of viable cells vs. tumour 
characteristic size x (^m), assuming X = 100 ^m. The set of viable 
cells becomes surface-like as the tumour grows. Therefore nearly all activity 
is eventually confined to a thin layer between the bulk of the tumour and 
the nourishing medium (when in vitro), or between bulk and blood 
vessels {in vivo). 

tumours, in spite of their recognized biological differences. These 
mechanisms are necessarily related to parameters X and A. 

Discussion 

We discuss the role of parameter first, starting from the remark that 
the volume taken by live cells is^ 



Va{x)=F{x)V{x)-- 



32 



-Ax^ 



(9) 



Eq. (9) shows that the volume of live cells scales as Va{x) ~ x^, with 
a size- dependent fractal dimension 



D-- 



dlnVa 
d\nx 



2>k + x 



(10) 



and that I determines this fractal dimension, which interpolates 
between D = 3 at small x, and D = 2 at large x (see Fig. 3). This 
behaviour is not associated to any irregularity of the external tumour 
boundary surface, but rather to the dispersion of live cells in the 
tumour tissue, as viable cells are mostly concentrated close to the 
tumour surface in tumour spheroids or to blood vessels in vascular- 
ized tumours, and gradually die away in the deeper layers \ This 
fractal behaviour may remind one of the nutrient distribution net- 
work of the West, Brown, and Enquist^^ (WBE) model, which has a 
fractal-like subdivision, down to a "minimum tube size", and leads to 
a scaling law for the basal metabolic rate of living organisms which is 
a power-law with exponent 3/4. However, in the case of solid 
tumours we do not find a single power-law, but rather a curve which 
interpolates between two different power-laws, Eqs. (5) and (6). 
Indeed, the concept of "minimum tube size" of the nutrient distri- 
bution network introduced in the WBE modeP^ does not apply here, 
since the data shown in Fig. 2 include both avascular and vascularized 
tumours: it is the very fine-grained fractal behaviour associated to the 
distribution of live cells - and therefore the value of parameter X - 
that actually determines the observed metabolic scaling in the case of 
tumour spheroids and solid tumours. 

We already noticed in reference 1 that 1 somehow summarizes the 
interplay of different microscopic processes - most notably, the dif- 
fusion of nutrients and toxic metabolites, and the biomechanical 
mixing of cells - and thus the exclusive dependence of the scaling 
law (8) on X seems to downplay the role of the other parameter A, 
which is related to the overall shape of the tumour. As remarked 
above, A is one of the factors that set the timescale of tumour growth 
and has an importance of its own, as it determines both the total 



tumour volume, V = Ax', and the tumour surface area, S ~ 3Ax^. It is 
also important to stress that S corresponds to the boundary between 
the bulk of the tumour mass and the non- cancerous environment, 
and this includes the interface between tumour and blood vessels, 
even where they penetrate the tumour mass (see Fig. 4). 

In the case of highly vascularized tumours we expect to find large 
values for A, and indeed the human tumours included in Fig. 2 yield 
values of A much larger than those found in avascular in vitro models 
(see Table 5 in Methods). 

The wide-range generality of scaling law (7) - or its "doubly nor- 
malized" version (8) - is noteworthy, and points to a common mech- 
anism that drives tumour growth and bypasses the details of 
vascularization and the pecuharities of different histotypes. More 
accurate and specific experiments could establish this law with 
greater precision, and eventually lead to a clinical scale of tumour 
progression from measurements of glucose consumption. 

Methods 

Here we provide detailed information on selection and analysis of the experimental 
data that we used to test the metabolic model for solid tumours. We discuss the 
reasons for preferring and choosing a limited set of data among the many available in 
the literature. Basically, this choice depends on how experimental data have been 
obtained and reported. Eventually, we select a list of 83 independent observations, a 
number that we think it is reasonably large to test model behaviour. 

Issues in the selection of experimental data. In spite of the large number of available 
data, their usability is constrained by the following requirements: 

• the scaling law (7) holds for individual nutrients as well as for the total nutrient 
consumption, and c depends on the specific nutrient addressed by the model as 
well as on the cell line. Therefore we concentrate on glucose, because glucose data 
are those that are more readily available. 

• The variability of tumour shape means that A is variable as well, and therefore a 
simple fit of expression 

3^+(y/A)'/^ 

cannot determine either X or A, but only their product XA. This means that we 
must fix the constant ciVc and this is only possible when the cell-specific para- 
meters are known: in other words, both cell- specific data and whole-tumour 
growth data must be available, and they must refer to the same cell line; 




Figure 4 | Schematic illustration of a vascularized tumour (not to scale). 

Live cells (LC) are distributed along the tumour vasculature (TV) itself and 
on the tumour boundary where it is close to blood vessels in the normal 
tissue (NT; blood vessels in normal tissue are not shown). The value of A is 
determined by the interface between the bulk of the tumour and the non- 
cancerous environment, and this includes the interface with tumour blood 
vessels. Areas of tumour quiescence/necrosis (N) are also shown. 
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Table 1 | Curated 9L spheroid data. The final line in the table 
reports the corresponding mean values, standard deviation and 
coefficient of variation of cell volume and glucose consumption rate 

n. Volume (lam^- 10^) fi (|imol-min-^ • 1 0"^) 



1 


0.153 




0.606 


2 


0.241 




0.556 


3 


0.183 




0.728 


4 


0.424 




1.122 


5 


0.544 




1.320 


6 


0.833 




1.586 


/ 


u.ooy 




1 .ouo 


8 


0.983 




3.721 


9 


1.152 




4.236 


10 


1.262 




5.548 


11 


1.441 




6.385 


12 


4.363 




6.263 


13 


5.290 




7.367 


14 


7.028 




12.781 


15 


8.768 




11.358 


16 


10.842 




21.715 


17 


1 1.739 




27.243 


18 


13.782 




27.259 


mean 


3.871-10^ M 


7.867- 


1 0"^ (lamo 


a 


4.580-10^ M 


8.885- 


1 0"^ (lamo 


CV 


1.183 




1.130 



• for the same reasons, we discard data expressed in arbitrary units. This rules out 
many data like those obtained with uncalibrated imaging techniques, those 
obtained with radioactive markers without any indication of the specific activity 
of labelled compounds, and data normalized per gram of tumour tissue without 
any reference to the total amount of tissue used. 

After an extensive search of the scientific literature, we found a limited set of data 
on glucose uptake in spheroids and solid tumours that fulfil these requirements. In the 
following we discuss the processing of these data to obtain the final dataset of 83 
independent data points used to test the scaling model. We express all data in the 
following units: [inr' for volumes; |imol/min for the rate of glucose uptake. 

While it is common to find published values of the rate of glucose uptake with its 
corresponding standard deviation in cell cultures, it is unusual to find more than one 
measurement of /i in solid tumors, and often there is no estimate of its standard 
deviation. Thus, instead of using the straightforward error propagation formula 

in most datasets we are limited to its more restricted version 

t7,l«^t7, (12) 



where we make the minimal assumption cr^ = 0. We also assume that the cell volume 
distribution is symmetric, and we take the median of equivalent to the average. 

9L cells and spheroids - data from ref. 5. 9L cells are from an established rat 
glioblastoma cell line. These cells form three-dimensional spheroids when grown 
under selective culture conditions and the data for both cultured cells and spheroids 
are available^; in particular we take the data from Figs. 1 and 4 in ref 5, respectively. 

The maximum reported glucose consumption by 9L cells in culture is c = 1.33 ± 
0.042 pg cell"^ min~^ and the median cell volume = 1200 [inr', thus rj = (6.15 ± 
0.194)-10"'^ lamol min"^ \im-\ 

The glucose concentration in spheroid culture media was measured by enzymatic 
methods, and the consumption rate was estimated by measuring glucose depletion at 
given time intervals. Care was taken to correct the glucose uptake due to single cells 
present in culture, and then data were normalized per spheroid, and these curated 
data are listed in Table 1. 

Ratl-Tl and MRl cells and spheroids - data from refs. 6-7. Ratl-Tl cells derive 
from T24Ha-ras-transfected rat embryo Rati fibroblasts, and MRl cells from myc/ 
T24Ha-ras-co-transfected rat embryo MRl fibroblasts. Both engineered embryo 
fibroblasts acquire a tumorigenic phenotype and the cells can form spheroids under 
appropriate culture conditions. Glucose consumption rates for these cell cultures 
have been measured and the data, already normalized per cell volume (i.e. rj in our 
notation), have been reported, see Table 2 in ref 6. The reason for data normalization 
stems from the careful experimental observation that the two cell types change their 
volumes as the consequence of the transfection with different oncogenes, and thus the 
normalization is required for proper quantitative comparison. The rates of glucose 
consumptions for exponentially growing Ratl-Tl and MRl cells are rj = (3.48 ± 
1.38) -10"'' lamol min"i [im'^ for Ratl-Tl cells and t] = (8.94 ± 0.3) -10"'^ |amol 
min"^ [im~^ for MRl cells. 

Mean glucose turnover in Ratl-Tl and MRl spheroids was measured photome- 
trically by quantifying the decrease in metabolite concentration per time interval in 
culture media. Values of fi and can be obtained from Fig. 4 in ref 7, where they are 
plotted as functions of spheroid diameter (in |im) and are expressed in units of (moles 
of glucose uptake) s'^cell"^* 10"^^, i.e., on a single-cell basis. Moreover in the same 
paper the authors also report careful measurements of the number of live cells per 
spheroid, and therefore the glucose rates per spheroid can be calculated as well - 
indeed, ref 7 is a beautiful example of accuracy in data presentation. We take the data 
from Fig. 2 in ref 7, and we fit them with the function /(x) = aix^^ , where x is the 
spheroid diameter. The results are shown in Fig. 5 and Fig. 6. We use the estimated 
values of ai and (X2 to calculate the number of live cells in spheroids of given diameter 
and the overall glucose uptake rate per spheroid. Finally, we assume a roughly 
spherical spheroid shape to compute spheroid volume; the curated data are listed in 
Tables 2 and 3. 

Human tumour cells and human primary and metastatic solid tumours - data 
from refs. 8-12. Glucose uptake rates in human solid tumours grown in vivo were 
measured by Kallinowski and collaborators^. In Fig. 1 of their paper, these authors 
report measurements carried out with a number of human primary and metastatic 
tumours xenografted into athymic, T-cell deficient rats. Glucose concentrations were 



10^ 




Spheroid diameter (//m) 

Figure 5 | Number of live Ratl-Tl cells per spheroid as the function of 
spheroid diameter. Symbols refer to experimental data from Fig. 2 in ref 7. 
The line is the best fit with a power-law function (see text). The fit returns 
the following estimates for parameters: ai =0.0074 and 0(2 = 2.485. 



Table 2 | Curated Rati -Tl spheroid data. The final line in the table 
reports the corresponding mean values, standard deviation and 
coefficient of variation of cell volume and glucose consumption rate 

n. Volume (jim^- 10^) fi (|imol-min-^ • 1 0"^) 



1 


0.020 


0.082 


2 


0.051 


0.162 


3 


0.054 


0.202 


4 


0.179 


0.478 


5 


0.337 


0.927 


6 


1.015 


2.703 


7 


2.175 


4.518 


8 


5.235 


9.338 


9 


14.715 


24.956 


mean 


2.642-10^ M 


4.818-10-^ (^imol 


a 


4.835-10^ M 


8.140-10-^ (^imol 


CV 


1.830 


1.690 
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Figure 6 | Number of live MRl cells per spheroid as the function of 
spheroid diameter. Symbols refer to experimental data from Fig. 2 in ref. 7. 
The line is the best fit with a power-law function (see text). The fit returns 
the following estimates for parameters: ai =0.00268 and 0(2 =2.592. 

measured enzymatically in the blood flow of perfused tumour tissues by inserting 
cannulas into both the supply artery and the drain vein, and the kinetic decline of the 
metabolite concentrations in the blood before and after the tumour tissues was 
carefully used to compute the rates of glucose uptake by tumours. 

Data are given in units of |j,mol of glucose per gram of wet tumour tissue per 
minutes and are plotted as the function of the tumour wet weight (in grams). Since the 
mass density of tumours is approximately 1 pg |am"^, we infer both glucose uptake 
rates and tumour volumes from these data. The curated data are listed in Table 4. 

The human tumours studied in ref. 8 included breast, uterine and ovarian carci- 
nomas, melanomas, thyroid carcinomas and lung carcinomas. Unfortunately, they 
did not investigate glucose uptake in isolated cells nor did they measure cell volumes. 
We have therefore mined the scientific literature for data obtained with established 
human carcinomas cell lines in vitro. As we shall see below, the measurements do not 
vary much between different cell lines and this allows us to estimate an average rj value 
for human cells. 

In Table 3 of their paper^, Mazurek et al. report the glucose consumption rates of 
two human breast carcinoma cell lines, MCF7 and MDA-MB-443; the values are: 
CMCF7 = (43.8 ± 0.5) nmol (10' cells)"' hour" ^ Cmda = (10.6 ± 0.9) nmol (10' cells)"' 
hour"'. The measured volumes of MCF7 (4000 fL = 4000 |am"') and MDA 
(approximately 7000 \im~^) cells can be found in refs. 10 and 11, respectively. Thus, 
we calculate rjMCF? = (1-83 ± 0.01) -10"'^ |imol min"' [Lm~^ and vjmda = (2.5 ± 
0.21) -10"'^ lamol min"' |am"l 

On the other hand, Aykin-Burns et al. studied glucose consumption with human 
colon carcinoma cells and normal non-immortalized colon cells'^. They find the 
following values for the three tumour cell lines HT29, SW480 and HCTl 16: Cht29 = 



Table 3 Curated MRl spheroid data. The final line in the table 
reports the corresponding mean values, standard deviation and 
coefficient of variation of cell volume and glucose consumption rate 

n. Volume (lam^- 10^) fi ()imol-min"^ • 1 0"^) 


1 


0.035 


0.260 


2 


0.086 


0.553 


3 


0.291 


1.218 


4 


0.869 


2.821 


5 


1.814 


5.089 


6 


5.521 


10.841 


7 


12.989 


20.981 


8 


23.908 


34.816 


9 


30.613 


39.933 


mean 

a 

CV 


8.458-10^ M 
1.155-10^ M 
1.366 


1.294-10"^ (^imol-min"^) 
1.538-10"^ (^imol-min"^) 
1.188 



Table 


4 Curated human tumour data. The final lines in the table 


report the corresponding mean 


values, standard deviation and 


coefficient of variation 


of cell volume an 


d glucose consumption rate 




Volume 


fi 




Volume 




n. 


(^im^-lO^^) (^mol-min"i) 


n. 


(^im^-lO^^j (^mol-min"^) 


1 


0.329 


0.707 


25 


1.354 


1.315 


2 


0.349 


0.360 


26 


1.480 


1.122 


3 


0.430 


1 .744 


27 


1.525 


1 .937 


4 


0.507 


1.167 


28 


1.525 


2.997 


5 


0.522 


1.111 


29 


1.334 


4.137 


6 


0.579 


1.116 


30 


1.796 


2.374 


7 


0.597 


1.01 1 


31 


2.023 


2.545 


8 


0.597 


0.972 


32 


2.147 


2.755 


9 


0.579 


0.789 


33 


2.212 


2.087 


10 


0.624 


0.91 1 


34 


2.766 


3.088 


1 1 


0.662 


0.967 


35 


2.935 


4.414 


12 


0.672 


1 .549 


36 


3.307 


3.278 


13 


0.703 


1.652 


37 


4.013 


6.865 


14 


0.693 


0.898 


38 


4.196 


6.064 


15 


0.693 


0.837 


39 


3.954 


4.639 


16 


0.769 


1.167 


40 


4.013 


3.675 


1 7 


0.906 


1.335 


41 


4.588 


4.415 


18 


0.990 


1.389 


42 


5.246 


3.979 


19 


1.051 


1.596 


43 


6.462 


5.917 


20 


1.067 


1.937 


44 


7.388 


5.716 


21 


1.149 


2.259 


45 


9.237 


8.210 


22 


1.202 


2.160 


46 


1 1 .377 


1 1 .504 


23 


1.099 


1.424 


47 


14.654 


9.673 


24 


0.933 


0.715 










Volume 




Glucose consumption 


mean 


2.494-10^2 




mean 


2.819 (jimol-min"^) 


a 


3.008-10^2 


(^m^) 


a 


2.470 (iimol- 


min ^) 


CV 


1.206 




CV 


0.876 





(9.0 ± 1.0) lamol (10' cells)"' day"', Csw48o = (17.0 ± 3.0) |amol (10*^ cells)"' day"' 
and CHCT116 = (21.0 ± 3.0) |imol (10' cells)"' day"'. 

To the best of our knowledge, measurements of cell volume for all these cell 
lines are not available. However, using sensitive fluorimetric methods. Tan et al. 
measured and reported volume values for a number of human carcinoma cell 
lines, among which SW480 cells and other colon carcinoma cells'''. In all cases cell 
volumes are close to 2 pL (i.e. 2000 |am"^), and taking this value we find 
r}HT29 = (3.12 ± 0.35) -10"'^ |amol min"' ^lm"^ rjsw48o = (5.9 ± 1.05) 
•10"'^ |imol min"' |im"^ and rjHcriie = (7.3 ± 1.05) -10"'^ [imol min"' [Lm~^. 

Finally, taking into account all the above values we find the average value 
rj = (3.68 ± 1.3) -10"'^ |amol min"' |am"l 



Table 5 | Fits of experimental data with Eq. (7): estimated para- 
meter values and goodness-of-fit statistics. Columns 1-3 are self- 
explanatory, column 4 is the chi-square per degree-of- freedom in 
the fit. V^hile the values of X are restricted to a rather narrow range, 
which is in line with the expectation of biological variability 
between different cell lines, the values of A cover a very wide 
range. This large variability is justified by the different geometries 
of the tumours in the list: while spheroids have shapes that never 
depart too much from the sphere - even when they have very 
irregular edges - real, vascularized tumours have complex geo- 
metries and the total area in contact with blood vessels is very 
large, hence the large average A value reported for human 
tumours 



Dataset 

9L spheroids 
Rat-Tl spheroids 
MRl spheroids 
Human tumors 



A (|im) 

93.93 ± 0.08 
123.56 ± 1.14 
90.95 ± 0.001 
99.40 ± 5.46 



2.75 : 
20.74 : 
0.76 : 
9648.2 : 



0.87 
2.27 
0.04 
: 0.02 



1.1 1 
1.54 
1.20 
1.02 
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Figure 7 | Nonlinear fits of the four data-sets with Eq. (7): symbols, 
experimental data; lines, nonlinear fit. The colours of both symbols and 
lines refer to different experimental data-sets and are as follows: blue, 9L 
spheroids; red, Ratl-Tl spheroids; yellow, MRl spheroids; green, human 
tumours. 



Estimates of k and A. We fit each selected experimental dataset with Eq. (7). The 
results are shown in Fig. 7. The corresponding values of parameters I and A, as well as 
the goodness-of-fit statistics, are given in Table 5. 

Data fitting. The data analysis is rather simple in this case, where data are 
comparatively few and the model is very smooth. To this end we used the 
NonlinearModelFit instruction which is built in Mathematica vs. 9.0 (Wolfram 
Research, Inc.). 
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